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Abstract 

We prove rigorous lower and upper bounds for the mass gap of the ferromagnetic spin 1/2 
XXZ chain. The two bounds coincide asymptotically in the Ising limit A — oo. Near the 
isotropic point, A = 1, the estimates are good enough to determine the critical behaviour of 
the mass gap unambiguously. The derivation does not rely on exact solutions. 
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1 Introduction 

The main purpose of this paper is to estabhsh rigorous upper and lower bounds on the 
spectral gap of the one-dimensional ferromagnetic spin 1/2 XXZ Heisenberg chain. See 



below in Theorem ^]T]for the explicit formulae. The upper bound is generally assumed to be 
the exact result. The lower bound has essentially the same behaviour and is reasonably close 
to the upper bound. It is good enough to determine unambiguously the critical exponent a 
of the model. We hope to make clear that our method relies only on certain properties of 
the ferromagnetic XXZ chain and not on the exact solution. We expect that the method will 
be useful in the study of other models that share the same general properties without being 
exactly soluble. In fact, similar ideas were already applied to a class of quantum spin chains 
with discrete symmetry breaking in ||l[] . We refer the reader to Section ^ for a discussion of 
the method and further references. 

The ferromagnetic regime of the XXZ chain has not been studied as extensively as the 
antiferromagnetic and the critical regimes. One could think that this is because it is less 
interesting or less challenging. We think this is only partly correct. In fact, the literature 
makes it very clear that the model is not so well understood as one might deduce from the 
fact that it is "Bethe Ansatz soluble". The more careful practitioners of the art do not 
neglect to point out that the validity of the Bethe Ansatz solution relies on the so-called 
string hypothesis, which remains unproven till now. In fact, it is known that the string 
hypothesis cannot universally hold, i.e., it is violated for some finite chains 0. 

Proofs of the completeness of the Bethe Ansatz based on combinatorial arguments (count- 
ing the number of solutions of the Bethe equations) always assume the string hypothesis . 
The proof of Yang and Yang of the validity of the Bethe Ansatz for the ground state 
of the XXZ chain only covers the range — oo < A < — 1, which is the complement of the 
regime studied in the present work. For the case A = 1 a form of completeness in the 
thermodynamic limit was shown in [Q. 

In view of the general uniqueness theorem for the Gibbs state of one-dimensional 
quantum spin models with translation invariant finite range interactions, there should be 
no surprises in the finite temperature behaviour of the XXZ chain. Yet, there are still 
unresolved questions about the low temperature behaviour of the specific heat of the model 
in the ferromagnetic region (See and p. Chapter 6] for a discussion.) 

These unresolved questions clearly demonstrate the need for rigorous arguments. The 
arguments presented in this paper do not address all of them, but they provide unambiguous 
information on the behaviour of the mass gap. 

The XXZ-Hamiltonian for a finite chain of L sites, including the special boundary con- 
ditions that we consider, is 

i/f = A{A){Sl - Sf) - E UslSl^, + SlS!^,) + SlSl_,, (1.1) 

x=l ^ 

where S^,a = 1, 2, 3, are the usual 2x2 spin matrices (with eigenvalues ±1/2) acting on the 
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site X, and A{A) = i^y 1 — 1/A^. The two Hamiltonians corresponding to the positive and 
negative choice of ^(A) are obviously unitarily equivalent by left-right symmetry. Unless 
explicitly mentioned we will always refer to the Hamiltonian with the positive choice for 
A{A). The boundary conditions and normalization of (|1 . 1|) are natural for the following 
reasons. First of all they make the ground state degeneracy equal to L + 1 for all A > 1. 
This property can be explained in terms of a quantum group symmetry that the Hamiltonian, 
with these particular boundary terms included, possesses . The normalization is such that 
one can consider the limit A — > cxo without difficulty. In this limit the model becomes the 
ferromagnetic Ising chain with a boundary term that allows for ground states with a kink, 
i.e., for any site x in the finite chain the configuration with all spins to the left of x up (|) 
and all spins to the right of x down (|), is a ground state. Obviously there are L — 1 of 
such kink states. Together with the two translation invariant configurations this yields L + 1 
ground states. As we will see in Section |^ the boundary terms also make the computation 
of the GNS Hamiltonians of the infinite chain immediate. 

In the isotropic limit (A = 1) the L + 1-fold degeneracy is the dimension of the spin 
L/2 representation of SU(2). Note that the boundary terms vanish for A = 1. In the 
thermodynamic limit {L oo) all translation invariant ground states are states of perfectly 
aligned spins. No non translation invariant ground states are known. 

If A > 1, there are four different classes of known ground states of the model on the 
infinite chain, which could be called up, down, kink, and antikink. They consist, respectively, 
of the state with all spins |, the state with all spins |, an infinite number of states in which 
the spins are | at — oo and | at -|-oo, and an infinite number of states in which the spins are 
I at — oo and t at -|-oo. The infinite degeneracy of the ground state in the latter two sectors 
corresponds to the possible choices for the location of the kink or antikink, i.e., the location 
on the chain where the spins turn over form up to down in the first and down to up in the 
second case. The kinks are strictly speaking located at a single bond only in the Ising limit 
(A oo). For 1 < A < oo the ground states are not described by a single configuration 
because of the quantum fiuctuations, but the kinks, respectively antikinks, are quasilocalized 
for all 1 < A < oo. We refer to Section ^ for more explicit properties of the ground states. 

Our main results are the upper and lower bounds on the gap in the spectrum of the 
Hamiltonian in the thermodynamic limit in the aforementioned four superselection sectors. 
Each of these sectors corresponds to a different representation of the observable algebra of 
the system. In these representations the Heisenberg dynamics of the model is generated by a 



densely defined self-adjoint, non-negative definite operator H. Theorem |1.1| below refers to 
the gap above zero in the spectrum of this operator. Alternatively, we can define the gap of 
the model with respect to an infinite volume ground state ( • ) as the largest constant 7 > 
such that for all local observables X 

{X*H'^X) = {X*[H\X]) > -f{X*[H^,X]) = -f{X*H^X) (1.2) 

By a local observable X we simply mean here a polynomial of a finite number of spin matrices. 
All local observables obtained from spin matrices with site index x in a given finite subset 
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A of the chain form an algebra which is denoted by ^a. The inequahty ( |1.2| ) expresses the 
property that the Hamiltonian is larger than 7 on its range, i.e., on the othogonal complement 
of the space of ground states. If the ground state is non-degenerate in the representation 
under consideration, the following equivalent inequality is more customary: 

{X*[H,X])>^{{X*X)-\{X)\') (1.3) 

In both (|1.2|) and ( |1.3|) the commutator has to be interpreted as the limit of successive 
commutators 

[H^, X] = hm ■ ■ ■ hm [Ha,, [i^A,.,, ■ ■ ■ [i^A,, X] ■ ■ ■]] (1.4) 

Ap— >oo Ai— ►oo 

which, due to the fact that the interaction is of finite range, is a local observable for all local 
X. 

Theorem 1.1 In each of the sectors described above as up, down, kink, and antikink, the 

infinite volume gap 7 satisfies 



for A > 1 and where 



riow(A) < 7 < rupp(A) (1.5) 

r„pp(A) = l-A-i (1.6) 

riow(A) = f 1 - A 1 - Ji-^ I (1.7) 



For A = 1, the difference •jl between the lowest and second-lowest eigenvalue of the Hamil- 
tonian Hl of a finite chain of length L, L > 2, satisfies 



^ <a-\l-T^r<lL<^n\j- + ±) (1. 



4L2 - ' V L - r - - L3' 

Two other parametrizations common in the literature are given by A = cosh A = (/x + 
fi~^)/2. The range 1 < A < +cxo corresponds to < A < -\-oo and 1 > /x > 0. The parameter 
fi is often denoted by q. The most common Hamiltonian is A times H^^^ defined in (|1 . 1| ) 



(up to the boundary terms). In terms of the parameter A the bounds of Theorem |1.1| are 

AFupp = coshA-1 (1.9) 



AFiow = cosh A (1 - Vl -tanhA ) (1.10) 



Near the isotropic ferromagnet (A = 1) the upper and lower bound both behave linearly, 
with a slope 2 and 1 respectively. This fixes the critical exponent a which governs the 
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behavior of the gap (as well as the low-temperature behaviour of the specific heat) near 
A = 1 to be equal to 1, which is in agreement with the exact result of Johnson, Krinsky and 
McCoy |10]. Near A = -|-cxo the lower bound behaves as 1 — \/2/A. 

Figure 1 gives an idea of the difference between the bounds Flow and Fupp. Fupp, here 
derived by a variational argument in Section HI is the exact solution given by the Bethe 



Ansatz pT] , |I2| |I3| . Also for the finite volume estimates at A = 1 , there are exact expressions 
for the coefficient of l/L^ g. 

The advantages of the approach in this paper are: 1) The method can be applied also 
to Hamiltonians that cannot be explicitly diagonalized; 2) As the the treatment is short, 
transparent, and completely rigorous, it should also deliver a better insight in the determining 
properties of the gap; 3) It is possible to obtain useful bounds down to the critical point 
(here A = 1), whereas this is usually not possible with other rigorous methods such as, e.g., 
the polymer expansion technique of fir 



2 The ground states of the XXZ chain 

Only these aspects of the ground states of the XXZ chain that have direct relevance to our 
estimates and understanding of the spectral gap of the model will be presented here. A more 
detailed analysis can be found in ||l6l, and various aspect of the ground states have been 



discussed in the literature (see e.g. and the references therein). It should be mentioned 
that a full analysis of the ground state problem for the infinite chain has not been achieved 
yet. Below we give a clear description of what is believed to be the complete set of ground 
states for the infinite chain. I am not aware of a rigorous proof that this is indeed the 
case. Loosely speaking one would obtain a description of the complete set of ground states 
by studying the thermodynamic limit with arbitrary boundary conditions. The difficulty is 
that a simple description of a sufficiently large class of boundary conditions is not available. 
Fortunately the statements in this article do not depend on the completeness of the set of 
ground states considered. 

For the study of the finite chains in this section we shall employ the special boundary 
conditions introduced in (|1 . 1|) . This choice of boundary conditions simplifies the study of 
the thermodynamic limit. It is also convenient to add a constant to the Hamiltonian to 
make the ground state energy vanish. This way, using the parameter /i, we arrive at the 
Hamiltonian, 

Hi = Hf""' + {L- l)/4 = E (2-1) 

x=l 

where h'^^^_^_i is the orthogonal projection on the vector 

e, = T^(/^in)-ur)) (2.2) 
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In terms of the spin matrices hi 2 is 

= -^-\SlSl + SfSi) - SfSl + ^ + A{A){Sl - S!) (2.3) 

with A defined following (|1 . 1|) . From the definition of it is obvious that h'^^x+i IT ' ' ' T) = 
for all a; = 1, . . . , L — 1. As is the sum of the h'^^^i, which are positive, this implies that 
the ground state energy of if^ is zero and that |t ■ • • T) is a ground state. For all < < 1, 
the ground state space (= ker H^^) is L + 1-dimensional. An explicit description of the L + 1 
ground states can be given in several ways. 

For all A > 1 the uniform states |T ■ ■ ■ T) and |i ■ ■ ■ i) are ground states of the XXZ 
chain. If A = 1 the L + 1-dimensional ground state space is the spin L/2 representation 



of SU(2). For all A > 1 and A = ^yl — l/A'^, the non- uniform ground states can be 
thought of as kink states, which are roughly described as the Ising kinks plus quantum 
fiuctuations. In this picture the degeneracy corresponds to the possible locations of the 



kink. For A = — lyl — l/A^ the kinks have to be replaced by antikinks, i.e., the roles of | 
and I spins have to be interchanged (or, equivalently, one can interchange left and right). 
We refer to [|l^ and 0] for more details and explicit expressions. 

In the thermodynamic limit the boundary terms disappear to infinity and the left-right 
symmetry of the model, broken by the particular boundary terms we have introduced, must 
be restored. It is therefore obvious that both the kink and antikink states appear as infinite 
volume ground states of the model. 

For our purposes the most convenient way to describe the space of ground states of a 
chain of length L is to introduce deformed raising and lowering operators which, together 
with the third component of the spin, generate the algebra (quantum group) of SU^(2). 
The concrete representation of SU^(2) is not left-right symmetric, and is different for the 
boundary terms that produce kink and antikink ground states. In fact the two mutually non- 
commuting representations of SU^(2) together generate the infinite-dimensional quantum 
affine symmetry algebra sl{2) that lies at the basis of the integrability of the model (see e.g. 
|T8|). We should stress, however, that a rigorous formulation of this infinite dimensional 



symmetry of the XXZ chain, has not yet been obtained. We will not use it here. 

In our computations we will not need anything beyond some basic facts of the representa- 
tion theory of SU^(2). We therefore restrict the discussion of the quantum group symmetry 
of the XXZ model to the bare minimum. One can think of the quantum group symmetry 
as a systematic way to construct operators that commute with the Hamiltonians H^. The 
parallellism with the usual arguments in the "theory of angular momentum" in quantum 
mechanics (representations of SU(2)) is so perfect that the reader will hardly notice the 
difference. 

For < /i < 1 define the 2x2 matrix t by 

t = (2.4) 
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and define as usual S"^ = ± iS"^. It is trivial to check that S"^ and t satisfy the following 
commutation relations 

tS^ = ^^^SH (2.5-a) 

[S+,S-] = ^^r^ = 2S^ (2.5-b) 
^ — /i 

They are just the SU(2) commutation relations in a disguised form. The remarkable fact 
is that there is a simple definition of the tensor product (coproduct of the quantum group 
or pseudogroup ||T9| , pO| ) of any two representations of the commutation relations ( p.5-a| )- 
( p.5-b| ), yielding a new representation. Here we only need the total-spin operators for a chain 
of L spins, which are given by 

L 

S[i,L] = Eli®---®'^'® ■■■II (2.6-a) 



x=l 



L 

^t,L] = ® Ix-+i®---Il (2.6-b) 



x=l 



L 

S\vM = Eli®---®'^^ ®^x+i®---^Z' (2.6-c) 



x=l 



where we used an index to identify the sites on which the tensor factors act. Note that, for 
L > 2, the operators depend on /i through t. One can easily check that the total "spin" 
operators as defined in ( p.6-a| )-( |2.6-c| ) commute with the interaction terms h^^^j^i and hence 

L 



with the Hamiltonian Hj itself 



3 Estimate of the mass gap for finite chains 



We begin this section with the derivation of a simple lower bound for the spectral gap of 
finite chains for Hamiltonians that share some of the basic properties observed in the XXZ 
chain (Theorem |3.2| ). A more general version of this estimate was given in where it was 
used to prove the existence of a spectral gap in arbitrary Generalized Valence Bond Solid 
chains with a finite number of ground states. As a strategy for obtaining lower bounds for 
the spectral gap of the generator of a spin dynamics, the method of proof is inspired by the 



work of Lu and Yau on the gap in the spectrum of the Glauber and Kawasaki dynamics 
of the Ising model. The ingredients that go into the estimate are not very different from the 
ones in 



22| and in fact similar elements underly the arguments in ||2^, ^ 



Here, we restrict ourselves to the simplest form of this estimate, which is sufficient for 
the application to the XXZ chain. 
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Consider an arbitrary spin chain of L sites and with Hilbert space TYl = ^x=ii^'^)x, 
where again we use the index x to associate the tensor factors with the sites in the chain. 
We assume that the Hamiltonian is of the following form: 



L-l 



Hl = X! ^x,x+i (3.1) 

x=l 

where h^^x+i is a translation of /ii^2, acting non-trivially only at the nearest neighbour pair 
{x, X + 1}. Assume furthermore that hi^2 > and that kei Hl ^ {0}. We will denote by 72 
the smallest nonzero eigenvalue of hi^2, i-e., the gap of H2. It is obvious that 

L-l 

ker ifi = Pi ker (3.2) 

x=l 

For an arbitrary subset A let G\ be the orthogonal projection onto 

ker h^,x+i (3.3) 

x,{x,x+l}cA 

For intervals [a,b], 1 < a < b < L, G[a,b] is the orthogonal projection onto the zero eigen- 
vectors of I]x=a hx,x+ii and Gs^^} = I for all x. From these definitions it immediately follows 
that the orthogonal projections Ga satisfy the following properties: 

G\2^Ai — Gj^^G\^ = G\^ if Ai c A2 (3.4-a) 

Ga.Ga, = Ga^Gai if Ai n A2 = (3.4-b) 

> 72(I-G[,,,+i]) (3.4-c) 

Define operators E^, 1 < n < L, on TY^ by 



Er, 



{ I - G[i,2] if n = 1 

G[i,„] - G[i,„+i] if2<n<L-l (3.5) 
G[i,L] if n = L 



One can then easily verify, using the properties (|3.4-a|) - (p.4-c|) , that {En \ 1 < n < L} is a 
family of mutually orthogonal projections summing up to I, i.e.: 

L 

= EnEm = Sm,nEn, ^ En = 1 (3.6) 

n=l 



The preceding paragraph applies directly to the XXZ chain. Next, we make a non-trivial 
assumption which we will verify for the XXZ chain later. 
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Assumption 3.1 There exists a constant e, < e < such that for all 1 < n < L — 1 

EnG 

[n,n+l] En < e^En (3.7) 

or, equivalently, 

\\G[n,n+l]En\\ < ^ (3.8) 



Note that, due to ( |3.4-a| ), G\n,n+i]En = G\n,n+i]G[i^n] — G[i^n+i\- This relates Assumption 



3TT| with Lemma 6.2 in ||2^, where an estimate for \\G\n,n+i]G[i^n\ — G\i^n+i] || is given for general 
Valence Bond Sohd chains with a unique infinite volume ground state. The same observation 
also implies that [G[„,n+i], Gii^^j] = [G[,„^„+i], which, if ( |3.8| ) holds, is bounded above in 
norm by 2e. 

The next theorem is a special case of Theorem 2.1 in Just like Theorem 6.4 in 
it provides a lower bound on the gap of the finite volume Hamiltonians, but it achieves this 
in a slightly more efficient way. We will repeat the proof for the particular case stated here, 
because it is simple, short, and instructive. 

Theorem 3.2 With the definitions of above and under Assumption \3.i[ the following esti- 
mate holds for all ip satisfying Gii^^ip = 0, i.e., ip that are orthogonal to the space of ground 
states of Hl : 

(V^|iJ^^)>72(l-v^£)l^f (3.9) 

i.e., the spectrum of Hl has a gap of at /easi 72(1 — -\/2^)^ above the lowest eigenvalue, which 
IS 0. 

proof: 

From the properties ( pl6|) of the £"„ and the assumption that i]?/' = 0, it immediately 
follows that 

llV^f = XillK^f (3.10) 

n=l 

One can estimate in terms of {ip \ hn^n+i4^) as follows. First insert and the 

resolution {Em}'- 

L-1 

WEni^f = I (I - G[„,„+i])E„^) + I E EmGyn,n+l]En^) (3.11) 

m=l 



Using ( |3.4-a|) and (|3.4-b|) one easily veryfies that Em commutes with Gn,n+i if either m < n—2 



OT m > n + 1. In these cases EmG[n,n+i]En = G[n,n+i]EmEn = 0, because the En form an 
orthogonal family. By this observation we obtain the following estimate. For any choice of 
constants Ci, C2 > 0: 

\\En^f = (V' I (I - Gln,n+l])En^) + ((^„-l + ^n)^ | G^^.n+ll^n^) 
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< -L(^|(i_G'f„,„+i])^) + |(^|E» (3.12) 



where we have apphed the inequahty 



for any c > 0, to both terms of ( 3.11 ). The first term in the right side of inequality ( 3.12| ) 



can be estimated with the interaction using ( p.4-c| ). The third term can be estimated with 
It then follows that 

(2 - Ci - -)||^„^f - C2||(K-1 + ^n)^ir < K,n+li^) 

C2 Ci72 

The term containing En-i is absent for n = 1, and En'ip = if = L. We now sum over n 
and use ( p.lO| ) to obtain 

(2 - ci - - - 2c2)||^l|' < ^(V' I Hl^) 

C2 Ci72 

Finally put ci = 1 — eV2 and C2 = s/V2 and one obtains the estimate ( p.9|) stated in the 
theorem. 

■ 

We now return to the XXZ chain and prove a lower bound on the gap of the finite volume 



Hamiltonians ( |1 . 1| ) as an application of Theorem 



Proposition 3.3 For a spin 1/2 chain of length L, the Hamiltonians 



L-l ^ . 



H'^'^^ — A{S\ S\) y^{-^('S'^'5'^+i + SlS^j^i) + S^S^j^i -} (3.13) 



with A = i^yl — 1/A^, A > 1, have an L + 1-fold degenerate ground state with eigenvalue 
and their next largest eigenvalue, 7^, satisfies 



7L > (1 - W^^)' = (1 - Vl-tanhA)2 = ri,„(A) (3.14) 

V + /i 

where rio„(A) is defined in ( \l.Ti) . // A = 1 one has the lower bound 
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proof: 

Due to the reflection (left-right) symmetry of the interaction it is clearly sufficient to consider 
one sign of A in ( p.l3| ), say A > 0. The proof consists in giving a constant e for which 
Assumption |3]^ holds. In fact, for the model under consideration one can simply compute 
the quantity The spaces on which the G[a,b] project are described explicitly in 

Section |[ Here, we consider all operators as acting on Hl- The quantity ||G[n,n+i]-En|| is 
equal to C„ defined for n > 1 by 

a= sup (3.16) 

0#V6H„+i \\lp\\ 

It follows from ( p.4-a| ) that H-EnV^H^ = — HG'ii^n+ijT/'Ip. Therefore, the requirement 
Enip = ip can be expressed as 

G[i,n]ip = ^, G[i^n+i]i^ = (3.17) 

As long asl<n<L — 1, C„ does not depend on L. 

The representation theory of SU^(2) is isomorphic to the one of SU(2) The ir- 

reducible representations can be labeled by the half-integers s = 0, 1/2, 1,3/2, . . ., and we 
denote them by D^'^\ We will use a subscript to indicate the set of sites on which a partic- 
ular representation acts. D'--^'^ is reahzed on 'H2S = i<S> C^)^** as the space of ground states 
of H2s- These elementary facts are sufficient to determine the vectors satisfying (|3.17|) . 



Indeed, G'[i,„]V' = ^ implies that e d[1^^^ ® ofj^l^ C Hn+i- As /^("/^^ ® D^^/^) 
^{(n-i)/2) ^ ^({n+i)/2)^ G[i,„+i]7/^ = implies that ip G D({"-i)/2). 

The ratio of norms in ( |3.16| ) is invariant under the action of SU^(2), and is therefore the 
same for all ^ G Z}(("~i)/2). One shows this just as one would for a group representation. For 
completeness we include a detailed argument here. G[„,„+i] = I — |^^) (^^| commutes with 
Sfi n+i] ■ Therefore, it is sufficient to consider ip that are eigenvectors of S^^ . Starting 
form one such eigenvector we can generate all others by application of •S'j'^^^^j. Assuming 
that S^^_^^^ip 7^ and using the fact that G[n^n+i] also commutes with S^^_^_^ we find 



(3.18) 



As {S^,n+i]y = %,n+i]Sii^n+i] (§66 ( p.6-a|) and ( p.6-b| )), (5f^^„+i])*5j„_^_i]^ is an eigenvec- 



tor of 5*^^ with the same eigenvalue as ip. This vector also belongs to the same irreducible 
representation and hence it must be proportional to ip. It follows that the ratio in ( p.l8| ) is 
the same as for ip. 

The computation of C„ is now straightforward, for we have to consider just one vector 
ijj E Z}(("^i)/2)^ e.g., one with S^^ = {n — l)/2. Any such vector is of the form 

n+l 

xlj = J2a,D, (3.19) 
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where denotes the ususal basis vector with all spins up except at the site x where the spin 
is down. Up to normalization the coefficients are uniquely determined by the conditions 
( p.lTp . A possible choice is 



^ = ^^ - J2 (3-20) 



A short computation shows that 



2 



(3.21) 



'^^"'"+^1^" = (/.-/.-T(/^ + /.-i) — ^^-^^^ 



and therefore 



y/i(l - /i2-2) 

a = / (3.23) 



which is increasing in n and C„ < l/-\/2 for all n > 1. We conclude that the assumptions of 



Theorem 3.2 are satisfied with 



e = < (3.24) 
\j fi + fi ^ 

This proves the theorem for < /i < 1. For fi = 1 one has to consider the limit 



e = limCL_i = ^Wl- (3.25) 



which is straightforward to compute. 



4 The infinite chain 

In order to be able to prove rigorous statements about the spectrum of the infinite chain 
we need to introduce the mathematical objects that define the infinite system. Although all 
interesting properties of the infinite chain can be expressed as results for limits of quantities 
defined for finite chains, the reverse is not true. Not all limits of finite chain quantities give 
interesting or even sensible statements about the infinite chain. By using a clean definition 
of the infinite system we will have no difficulty in sorting out the relevant statements about 
the thermodynamic limit of the XXZ chain. 
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Let the symbols TT)Ti)iT)ii denote the four superselection sectors of the infinite XXZ 
chain with A > 1, corresponding to up, kinks, antikinks, and down respectively. We can 
describe the GNS Hilbert spaces of these four superselection sectors as the so-called 



incomplete tensor products ||28| Ha/s, for a and /3 =| or |, defined by 



= U ® ^c.(3{y) 1 (4.1) 

A \x&A yeA= 



^<^M = \ m "n (4-2) 



where 

O J.A ^ / I") 

if y > 

We also define the vectors Qa/s as the infinite product vectors 

^af3 = <S> ^o.p{y) G (4.3) 

j/gZ 

Let A\ denote the local observables acting nontrivially only on the sites in the finite set A. 
Local observables X G act on Tiap in the obvious way, e.g., the spin matrices at the site 
X act on the x*'^ factor of the tensor product ( ^.1| ). From the definitions above it is clear 
that vectors of the form 

ilj = X^^p, XeU^A (4.4) 

A 

form a dense subspace of Tiap- Note that if a 7^ /3, VL^p is not the GNS vector representing 
one of the kink (or antikink) ground states. 

The mass gap of the infinite chain is a property which is defined with respect to a partic- 
ular ground state of the infinite system or, more precisely, with respect to a superselection 
sector. The Hamiltonian is represented on Ti^^ as the generator Hap of the Heisenberg dy- 
namics of observables acting on Tiap. The dense subspace of the vectors ip defined in ( ^^ ) 
is in the domain of -ffa/j, and the selfadjoint operator H^p is uniquely determined by the 
requirement 

H^pXn = lim [H^''^, X]^l^p (4.5) 

A-*Z 

We remark that Hap does not depend on boundary terms such as ^(5*^ — 5*^) added to the 
XXZ Hamiltonian for finite chains. It is well-known that Hap is a positive operator in 



general and in the present case this could not be more clear. An explicit formula for Hf^p is 

HapXQaP = X! K:'^x+l^^aP (4.6) 

{x,x+l}nA7^0 

where can be taken to be /i^^^^+i if (^P =1?! ii? oi' Ti- If C(P =iT ths sign of the 

boundary term has to be reversed. This is equivalent to replacing fi by fi~^. 
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The mass gap 'jafs is then just the gap above in the spectrum of Ha/s- A formula for 

7«/3 is 

7"/3 = _inf„ ^, , , I (4.7) 



As {k.eiHajs)'^ = ran Hai3 the infimum in (|4.71) can be taken over vectors of the form Hajjip, 
and it suffices to take ip of the form ( [4. 4] ) because they are a core for -ffa/j. 

There is no a priori reason why the spectrum of Ha/s should be independent of the 
superselection sector, i.e. independent of a/?. We already know that the multiplicity of the 
lowest eigenvalue is different: it is 1 for H-^^ and and infinite for H-^i and Hi^. Therefore, 
a priori, we should not expect 7^/3 to be independent of a/5. One can easily convince oneself, 
however, that 7q,^ = 7/3Q, and that 7|| = 7||. From a simple argument given in Section it 
follows that 7q,/3 < 7i-|. The upper and lower bounds that we will derive here are independent 
of ap. 

4.1 Proof of the lower bound of Theorem 



In order to prove the lower bound of ( |1.5|) we simply have to show that the lower bound 
on the finite volume gap obtained in Section ^ remains valid in the thermodynamic limit, 
irrespective of the particular zero energy ground state that we are considering. It is important 
that the finite volume gap estimates were obtained for the "correct" boundary conditions 
of (|1.1| ). More explicitly we show that for any choice of ajS and all local observabels X the 
following inequality holds: 

I X* H^f^X^la/s) > riow(A)(fia/3 | (4.8) 

If a = P, flaf3 is the vector representing the unique ground state of ifa/j. If a 7^ /?, Qa/s is 
not a ground state itself (except in the Ising limit A = 00). But all the kink (if a/? =ti) 
or antikink states (if af3 =it) are represented as vectors in the Hilbertspace Tiais defined by 
([4.1|), and together these vectors span kei Hap. A proof of these statements can be found in 

The inequality (|4.8|) follows from Proposition |3.3| when one observes that for X G 



{n^p I x*HlpXn^(,) = {n^fs \ x*iH^^,Yxn^p) (4.9) 

Obviously, X*{H'^_^^)^X G Aa±3- Therefore the expectation value in the right side of (^^ ) 
can be computed in the density matrix pa±3 which describes the state flap in the finite 
volume A ± 3. The same is true for the right side of ( |4.^ ). We conclude that it is sufficient 
to ascertain that 

Tr pA±3X*(i7j^±3)3x > rio„(A) TTpA±,X*{H^^,fX (4.10) 
which immediately follows from the finite volume gap estimate of Proposition piB. 
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4.2 Proof of the upper bound of Theorem 

First we argue that it suffices to prove the upper bound for It is obvious that the gap of 
will satisfy the same bound. For the gap of the model in the kink and antikink sectors 
we have an inequality which can be derived as follows. The translation invariant ground 
states can be obtained as weak limits of the kink or antikink states by letting the position 
of the kink (or antikink) tend to ±00. We then have 



< inf hm N"""P I '"v- \-A±.ij -y-qp/ ^2) 

- A,xeA.n^±^{n^^\r^(x*iH^^^,yx)n^p) ' 

A,xeAA (fin I X*{H^^syxn^^) ^ ' ' 

where r„ denotes the translation over n lattice units in the chain. It follows that 70/3 < 7||. 

For the proof of the upper bound it is convenient to present the dense subspace of ran Hai3 
formed by the vectors of the form (|4.4|) in a slightly different way. Observe that the spaces 
ker if^ C Tiap are decreasing in A. Therefore, in order to assure that a certain ip belongs to 
ran H^iB, it suffices to check that ip _L kerify( for some suitable A. 

We fix an interval [1, n] and introduce the usual spin wave operators X^, k = 27rm/n, m = 

0. . . . , n — 1, given by 

1 " 

X, = —Y^e''^S; (4.14) 

x=l 

The normalization and the allowed values for k are chosen such that 

{Q^^\X;XkQ^^)=Sk,i (4.15) 

The vectors ip we need for the upper bound are linear combinations of two spin waves, 

1. e. ip = {ciXk-^ + C2Xk2)^-\i- Due to ( [4.15| ) we have HV^P = |cip + |c2p. For any pair of 
distinct ki,k2, the coefficients Ci,C2 can be chosen such that G[i,n]V^ = 0, i.e. ip -L ^erifj^^j. 
This follows from the fact that ker Hj^^ contains exactly one vector for each eigenvalue of 
Sfin]- AH vectors Xkfl^i have S'pj^^j = (n — 2)/2. It follows that any two-dimensional space 
of vectors ip with fixed, distinct ki,k2 and arbitrary Ci,C2 must contain a ray _L kerifrj i. 



Hence the upper bound (|1.6|) of Theorem |0| can be proved by showing that 

inf =up<4M = r„„(A).l-i (4.16) 



which we do next. 
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From the definition ( [4.14| ) of the Xk it is clear that the only matrix elements of H^^ we 
need are the T^^y, 1 < x,y < n, defined by 

Tx,y = {^u I 5'+if[i5,n+i]'S'y ^Tt) = ^ {2A5^,y - So,,y-i - Sr,^y+i} (4-17) 



It is then easily seen that the sup^^^^ the left side of ( ^4.16| ) yields the norm of the 2x2 
matrix M(n, ki, /C2) with matrix elements 

M{n,h,k2h,j = M^{ki,kj) (4.18) 

where Mn{k, I), for k, I of the form 2Tim/n, is the function 

1 " 

M„(fc,/) = -Y. e-''^Tx,ye''y (4.19) 
^ x,y=l 

= 6k,i{l-A~^cosk) + {e'^ + e-'^)/{2An) (4.20) 

It is then obvious that min^kiM ^2)!! = 1 — A^"^. 

A similar calculation yields the L-dependent upper bound (^) in the case A = 1. 
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